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3/ 
In  1855  Fick  (4)~  put  diffusion  on  a  theoretical  basis  by 

showing  that  the  diffusion  of  substances  through  solutions  obeyed  the 

same  mathematical  laws  that  Fourier  had  derived  for  the  conduction  of 

heat  in  solids.   Since  the  time  of  Fick,  theoretical  and  experimental 

aspects  of  diffusion  have  undergone  spectacular  advances  both  in  single- 

and  multiphase  systems  (1,  9) 0   For  a  multiphase  system  such  as  a 

porous  material,  diffusion  takes  place  almost  exclusively  through  the 

fluid  phase  (gas  or  liquid  or  both).   The  solid  matrix  complicates 

matters  by  altering  the  diffusion  path  length  and  the  cross-sectional 

area  available  for  diffusion., 

The  diffusion  path  length  and  the  effective  cross  section  are 
extremely  complicated  properties  of  each  porous  material,  and,  in 
certain  cases,  both  may  be  functions  of  space  and  time,,   Accounting 
for  these  factors  in  diffusion  equations  is  usually  accomplished  by 
inclusion  of  the  two  parameters,  porosity  and  tortuosity,  which  are 
assumed  to  be  constants,,   The  manner  in  which  the  two  parameters  have 
been  used  in  diffusion  equations  has  not  been  universal.   Furthermore, 
the  adequacy  of  the  two  parameter  treatment  has  been  questioned 
(2.    6.  10). 

It  is  the  purpose  of  this  report  to  bring  together  the  several 
approaches  that  have  been  used  to  account  for  the  complications  caused 
by  the  porous  material  and  to  show  the  relationships,  and  some 
limitations,  of  the  several  concepts. 

The  discussion  is  confined  to  diffusion  equations  applied  to 
porous  materials o   Adequate  reviews  of  diffusion  processes  are  available 
(1.  8,  9) o   We  assume  isothermal  conditions  and  that  the  volume  of  the 
porous  medium  remains  constant. 
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FIRST  ORDER  EQUATIONS 

Quantitative  measurements  of  the  rate  at  which  a  diffusion 
process  occurs  are  usually  expressed  in  terms  of  a  diffusion  coefficient. 
For  brevity,  we  will  confine  our  discussion  to  one  dimension,  and  define 
the  diffusion  coefficient  as  the  rate  of  transfer  of  the  diffusing 
substance  across  a  unit  cross-sectional  area  divided  by  the  space 
gradient  of  concentration  of  the  substance  at  the  section.   Thus,  if 
q  =  Q/t  is  the  amount  of  substance  transferred  in  time  t,  A  the  cross- 
sectional  area,  c  the  concentration,  and  x  the  space  coordinate,  then 

q=-DA§f  [1] 

is  a  definition  of  the  diffusion  coefficient  D.   This  definition  can  be 
unambiguous  only  if  the  units  of  q,  A,  c,  and  x  are  carefully  specified 
(_8) .   As  Hartley  and  Crank  (_8)  have  noted,  concentrations  can  be 
expressed  in  a  number  of  different  ways,  thus  necessitating  that  all 
components  of  [l]  carry  an  index  to  indicate  the  frame  of  reference  used. 
In  general,  the  concentration  c  should  have  the  same  quantity  reference 
as  q,  and  the  volume  should  be  of  the  same  length  reference  as  A  and  x. 
The  use  of  different  frames  of  reference  for  components  of  the  same 
equation  requires  the  use  of  appropriate  conversion  factors.   This  is 
particularly  true  in  porous  materials  where  the  frame  of  reference  may 
be  the  entire  bulk  medium  (designated  by  the  subscript  m)  or  the  fluid 
phase  (designated  by  the  subscript  f ) .   In  some  cases  the  fluid  phase 
must  be  further  designated  as  gas  or  liquid  or  both. 

Writing  equation  [l]  in  reference  to  the  entire  porous  medium, 
we  have 

dc 

q  =  -  D  A  ^r-S    ,  [2] 

m  m  dx    ' 

m 

where  q  refers  to  the  quantity  (frequently  the  mass)  of  diffusing 
substance  per  unit  time  and  does  not  need  to  be  referred  to  the  medium 
or  fluid  reference.   If  two  or  more  substances  are  diffusing  simul- 
taneously, then  the  equation  can  be  written  for  each  substance  and  q 
would  carry  an  appropriate  index  (8)  . 

The  use  of  cm,  the  quantity  of  diffusing  substance  per  unit 
volume  of  porous  medium,  departs  from  the  more  commonly  used  definition 
of  concentrations,  but  it  is  the  only  correct  definition  that  can  be 
used  for  a  two-component  system  when  one  component  is  the  porous  medium. 
Its  use  is  not  without  precedent  (5,  13). 

In  many  systems  it  is  more  convenient  to  measure  the  concentration 

with  reference  to  the  fluid  phase.   If  e  is  the  relative  volume  (Vr/V  ) 

r  f  m 


occupied  by  the  particular  fluid  in  question,  then 

c   =  €C^      (0  <  £  <  1).  [3] 

m     r 

In  the  case  where  a  single  fluid  completely  occupies  the  pores  of  the 
medium,  e  is  the  porosity. 

The  correct  value  of  e  may  not  be  evident  in  every  experiment. 

For  example,  if  dead-end  pores  are  present  in  the  medium  (6 ,  7)  and  a 

steady  state  exists,  not  all  pores  contribute  to  the  diffusion  path. 

Hence,  the  total  volume  of  fluid  (V  )  may  not  be  the  correct  parameter 

to  use  in  calculating  e.   The  magnitude  of  the  error  involved  here  is  not 

known.   The  pertinent  literature  has  appeared  only  recently,  and  the 

data  are  largely  for  sandstone  models.   For  our  purpose  it  will  suffice 

to  mention  the  possible  complications  and  assume  that  e  is  uniquely 

given  by  V../V  and  is  a  constant, 
f  m 

Using  [3]  in  [2]  we  have 

i  "  "  V-  ST  [4] 

m 

where  D  =  D  e, 
p    m 

For  gaseous  diffusion  in  porous  materials,  Penman  (11)  defined 
the  diffusion  coefficient  as  given  by  equation  [4] .   This  convention 
has  been  widely,  but  not  exclusively  (14) ,  used  in  soils  literature. 

The  usual  definition  of  tortuosity  is  the  ratio  of  the  medium 
(or  externally  measured)  length  to  the  diffusion  path  length,  i.  e. 

xm  =  pxf     (  0  <  3  <  1).  [5] 

The  relationship  between  A  and  A,,  has  been  reported  in  several 
, .  rr  ,,  m      f 

different  forms: 

(a)  For  diffusion  paths  of  constant  cross  section  that  pass 
directly  along  a  line  perpendicular  to  the  surface  A  ,  Af  =  eA  . 
Several  investigators  have  used  this  form  (3,  11) . 

(b)  For  diffusion  paths  of  constant  cross  section  that  follow 

a  tortuous  path  through  the  medium,  the  following  argument  can  be  used. 
Since  e  =  Vf/Vm  =  Afxf/AmXm,  and  using  [5],  we  find  A^   =  Af(ge)-1,  a 
result  obtained  by  several  investigators  (2 ,  12 ,  15) . 

(c)  Diffusion  paths  that  have  a  nonuniform  cross  section  and 
also  follow  a  tortuous  path  through  the  medium  have  been  discussed  by 
Currie  (2).   He  argued  that  the  effective  area  available  for  diffusion 
may  be  reduced  by  constrictions  along  the  diffusion  path  and  treated  the 


simplest  form  of  this  case  by  assuming  that  the  area  along  the  diffusion 
paths  varies  sinusoidally. 

The  arguments  used  individually  in  the  above  three  cases  can  be 
combined  by  reverting  to  the  definition  of  the  diffusion  coefficient 
given  by  equation  [l].   Equation  [l]  was  written  for  an  infinitely  thin 
cross  section  of  area  A.   Applying  the  equation  to  a  porous  material,  we 
have  an  infinitely  thin  section  of  area  Am  in  which  the  area  available 
for  diffusion  is  eAm.   Consider  an  infinite  number  of  sections  so 
arranged  as  to  form  a  volume  A^m.   We  require  each  and  all  sections 
to  have  the  same  porosity  (e).it'   We  do  not  require  the  "holes"  in 
each  section  to  be  of  uniform  size  within  a  section  or  among  sections, 
nor  for  all  "holes"  to  have  the  same  space  coordinates  among  sections. 
We  require  that  the  "holes"  interconnect  to  form  a  continuous,  though 
tortuous,  path  of  varying  cross  section  through  the  thickness  xm.   In 
addition,  we  have  assumed  that  no  dead-end  pore  volume  exists  in  the 
medium. 

Consider  now  a  pore  shaped  so  that  its  area  increases  in  the  x 
direction,,   Since  each  infinitely  thin  section  must  have  a  porosity  e, 
another  pore  must  exist  whose  area  decreases  in  the  x  direction.   The 
net  area  of  the  innumerable  pores  in  the  thickness  Ax  would  be  constant. 

When  the  area  of  an  individual  pore  is  not  constant,  the  diffusion 
path  length  becomes  nebulous.   Within  an  individual  pore,  it  is  necessary 
to  take  some  average  value  for  the  path  length.   Thus,  instead  of  a 
tortuosity  factor  and  another  factor  accounting  for  the  nonuniformity 
of  area  in  individual  pores,  we  will  consolidate  these  factors  into  one 
"complexity"  factor  for  a  porous  material  and  designate  this  with  the 
symbol  a.   The  term  a   includes  (3  and  the  factor  accounting  for 
nonuniform  cross  section  in  an  unspecified  manner  and  is  associated 
with  length  in  the  same  manner  as  p.   We  must  associate  a  with  length 
because,  as  Ax  —>  0  the  effective  area  approaches  eAm,  and  the  diffusion 
paths  do  not  have  sufficient  length  to  exhibit  a  tortuosity  nor  a 
nonconstant  cross  section.   For  the  purpose  of  this  discussion,  we  will 
assume  that  a  is  a  constant.   Then 

x  =  a   x,,     (0  <  a  <  1)  [6] 

m      I 


and 


A  =  kc    (ck)"1   ,  [7] 


by  the  same  argument  used  previously  in  case  (b). 


4/   It  is  necessary  to  make  this  restriction  in  order  that  the 
area  term  in  the  diffusion  equations  can  be  considered  constant.   If 
it  is  not  constant  the  functional  relationship  between  eA^   and  xm 
must  be  specified  before  solutions  of  the  equations  can  be  obtained. 


Defining  Df  as  the  diffusion  coefficient  for  a  diffusing  substance 
in  bulk  fluid  (departing  from  the  more  conventional  notation  D0)  and 
using  [6]  and  [7]  in  [4],  we  find  the  relations  between  the  diffusion 
coefficients  to  be 

2 

D  =  D£  ea   ,  [8a] 

P    f 

and 

D  =  d.  a2       ,  [8b] 

mi 

In  applying  equation  [l]  to  the  steady  state  diffusion  of  a 
substance  in  a  porous  material  of  length  L,  we  need  only  to  require  that 
the  flux  q  be  constant  and  have  the  same  value  at  all  cross  sections 
within  the  system.   If  we  require  here,  as  we  did  in  the  discussion  of 
the  "complexity"  factor  above,  that  the  area  available  for  diffusion  be 
of  uniform  value  for  all  sections,  then  the  diffusion  coefficient  can  be 
evaluated  from 


qM„  -  -  d_  <%r    .  [si 


m      o  ox 
m 


if  the  concentration  is  known  as  a  function  of  distance,,   If  D  is 
constant,  then  [9]  reduces  to  the  familiar  form  p 

q/Am  =  D   (c    -  c  )/L  [10] 

m    p    r 1    tz 

where  Cf\   and  Cf2  are  the  concentrations  (with  reference  to  the  fluid) 
at  the  ends  x  =  0  and  x  =  L,  respectively.  Obviously  cm  and  Dm  could 
be  used  in  [9]  and  [ 10] .   The  experimental  conditions  usually  dictate 
the  frame  of  reference  to  be  used. 


SECOND  ORDER  EQUATIONS 

Consider  two  sections  a  distance  Ax  apart.   The  flux  entering  one 
section  is  q  and  that  entering  the  other  is  q  -  (dq/dx)  Ax.   Since  the 
rate  of  accumulation  of  substance  in  the  layer  AAx  is  A(dc/dt)  Ax,  we 
have 

If  the  cross  section  A  is  constant,  [ll]  becomes  the  familiar  second 
order  concentration-dependent  diffusion  equation,  viz 


If  D  is  also  constant,  [12]  takes  the  simpler  form 

ox 

Equation  [13]  is  known  as  Pick's  second  law.   The  first  order  diffusion 
equation,  usually  written  in  the  form  of  [l]  except  with  D  a  constant, 
is  known  as  Pick's  first  law.   This  terminology  is  probably  based  on  the 
fact  that  equation  [13]  is  a  second  order  equation  and  equation  [I]  is  a 
first  order  equation.   Fick's  original  equation  had  the  form  of  [11] 
with  D  constant  and  the  cross  section  as  a  function  of  distance  (4) . 
In  his  first  experiment,  he  utilized  equation  [13].   For  his  subsequent 
steady  state  experiments,  he  set  oc/ot  =  0  and  integrated  the  right- 
hand  side  to  obtain  c  =  ax  +  b.   This  is  of  the  same  form  as  [10]. 

If,  in  the  application  of  the  second  order  equations  to  porous 
materials,  the  cross  section  available  for  diffusion  varies  with  distance, 
then  equation  [11]  is  applicable.   The  use  of  [ll]  would  require  us  to 
know  exactly  the  functional  relationship  between  A  and  x.   For  the 
majority  of  applications,  it  is  again  convenient  to  assume  a  constant 
effective  cross  section  throughout  the  material. 

For  the  following  discussion,  either  equation  [12]  or  [13] 
would  serve  equally  well.   We  will  use  [13]  for  convenience. 


Equation  [13],  written  in  reference  to  the  entire  porous  medium, 


is 


oc       d  c 

"=  D  — =S  .  [14] 


ox 


m 


With  concentrations  expressed  on  basis  of  the  fluid, 


OCj.        O  Cj. 

1  =  D  — =i    .  [15] 


oT 


ox" 
m 


Thus,  regardless  of  the  concentration  reference,  the  diffusion  coefficient 

is  D  .   If  it  is  more  convenient  to  use  D  instead  of  D  ,  then 
m  p  nr 


dc£   D  d2c£ 

m 

Upon  integrating  equations  [14],  [15],  or  [16],  c  should  be 
expressed  in  both  integral  and  boundary  values  on  the  same  consistent 
basis.   If  the  integration  yields  q/A,  the  same  length  reference,  or 
appropriate  conversion  factors  for  c  and  A,  should  be  used. 

Comparison  of  equations  [2],  [4],  [15],  and  [16]  shows  that, 
without  properly  indexing  or  otherwise  identifying  the  diffusion 
coefficients,  confusion  can  arise  when  data  from  different  sources 
are  compared,  since  one  may  report  D  and  another  D  without  clearly 
pointing  out  which  was  used.   An  additional  complication  arises  if 
dead-end  pores  are  present  in  sufficient  numbers  to  cause  the  values  Djq 
(or  Dp)  measured  at  steady  state  to  be  different  from  values  of  the  same 
coefficient  measured  in  the  transient  state  (6,  7) . 

By  use  of  equation  [7]  in  [15]  and  [16]  the  same  relationship  is 
obtained  between  the  diffusion  coefficients  as  is  given  in  equations  [8a] 
and  [8b].   These  relationships  are  valid  only  when  a   and  e  are  constants. 
It  has  been  clearly  shown  by  Currie  (2)    that  Q^  (i.e.  Dm/Df,  his  figure 
2)  and  e  are  not  simply  related.   The  relationship  is  influenced,  among 
other  things,  by  particle  and  pore  shapes. 

Estimations  of  the  complexity  of  porous  material  have  been  the 
subject  of  numerous  investigations.   The  diffusion  coefficient  (either 
Djjj  or  Dp)  is  measured  and  related  to  Df  by  equations  of  the  form  [8]. 
For  the  ratio  Dm  (or  Dp)/Df  to  yield  a  significant  value,  Df  must  be 
known  with  an  error  less  than  or  equal  to  that  of  Dm  (or  Dp).   Values  of 
Df  taken  from  the  literature  may  not  always  be  an  appropriate  basis  for 
comparison.   For  example,  consider  the  diffusion  of  an  ionic  substance 
through  a  liquid  solution  held  by  a  porous  material  such  as  soil,  and 
the  diffusion  in  dilute,  unconfined,  single-salt  solutions.   It  is 
obvious  that  a  diffusing  material  encounters  radically  different  chemical 
and  physical  environments  in  the  two  conditions.   Thus,  chemical  factors, 
in  addition  to  the  physical  factors  enumerated  previously,  can  enter 
equations  [8a]  and  [8b]  (12) .   The  "complexity"  factor  is  then  even 
more  complex  than  that  indicated  in  the  earlier  discussion. 

It  is  doubtful  that,  for  diffusion  of  many  substances  in  soils, 
the  various  chemical  and  physical  factors  can  be  separated,  evaluated, 
or  enumerated.   The  data  reported  in  the  form  of  Dm  or  Dp  would  seem 
preferable  since  fewer  assumptions  must  be  made,  and  the  physical 
significance  of  D^  and  Dp  is  less  ambiguous  than  a  list  of  unseparable 
chemical  and  physical  factors. 

The  terms  Dm  and  Dp  are  equally  meaningful,  and  the  choice  of  one 
over  the  other  is  dictated  by  personal  preference  and  experimental  con- 
ditions,  whatever  frame  of  reference  is  used  to  denote  the  various  param- 
eters in  the  diffusion  equations,  it  should  be  described  unambiguously. 
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